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INTRODUCTION 


Oim  of  tho  central  goals  in  explosives  science  is  to  extrapolate  the  results  of  simple 
small-scale  testing  to  predict  the  outcome  of  complex  problems  of  practical  interest.  A 
remarkable  procedure  has  been  developed  by  C.  Forest  (Reference  1)  that  has  enabled  engi¬ 
neering  calculations  for  a  variety  of  problems  involving  shock-initiated  reaction  of  hetero¬ 
geneous  solid  explosives.  Forest  assumed  a  global  reaction  rate  model  for  the  decomposition  of 
the  explosive  ami  derived  a  method  to  calibrate  the  model  using  the  results  of  wedge  tests. 
The  result  is  popularly  known  as  the  "Forest  Fire”  burn-rate  model.  The  reaction  model  is 
then  used  in  one-  and  two-dimensional  hydrodynamics  computer  codes  where  the  details  of 
shock  wave  interactions  with  reacting  explosive  materials  can  be  calculated.  This  procedure 
has  been  successfully  applied  to  a  variety  of  very  practical  problems.  (See  References  2, 3,  and 
4,  for  examples.) 

In  the  selection  of  the  reaction  model  and  the  derivation  of  its  calibration  procedure, 
Forest  makes  an  assumption  about  the  processes  that  occur  in  the  explosive  decomposition 
that  is  contrary  to  the  widely  held  view  of  hot  spot  initiation  of  heterogeneous  explosives.  He 
also  makes  several  assumptions  about  the  flow  field  obtained  in  the  wedge  test  that  affect  the 
calibration  of  the  reaction  model.  These  assumptions  are  inconsistent  with  calculations  using 
the  burn  modal.  That  is,  when  the  Forest  Fire  burn  model  is  used  to  calculate  wedge  test 
results,  one  does  not  obtain  the  flow  field  assumed  in  the  calibration  of  the  model. 

In  spite  of  inconsistencies  in  the  theory,  one  is  faced  with  a  rather  impressive  history  of 
success  in  correlating  experiment  with  calculation  for  a  wide  variety  of  applications.  There¬ 
fore,  an  attempt  is  made  here  to  evaluate  the  effect  of  the  various  assumptions  of  the  Forest 
Fite  theory.  The  method  used  is  to  evaluate  the  reaction  rates  based  on  extreme  cases  of  the 
assumptions  and  to  evaluate  their  effect  on  results  of  numerical  calculations  of  the  wedge 
test 


FLOW  FIELD  ANALYSIS 


As  a  model  of  the  shock-to-deto  nation  transition  observed  in  the  wedge  test,  Forest 
proposed  a  hydrodynamic  discontinuity  within  which  some  fraction  of  the  solid  reactants  is 
converted  (burned)  to  gaseous  products.  In  one  extreme  case,  no  solid  is  burned  and  the 
discontinuity  is  the  familiar  shock  wave.  In  the  other  extreme,  all  of  the  solid  is  consumed 
and  one  has  a  detonation  wave.  In  the  following,  the  hydrodynamic  discontinuity  will  also  be 
called  a  reacting  shock  wave.  A  wave  diagram  of  an  accelerating  shock  wave  is  shown  in 
Figure  1.  At  a  distance  X*  the  wave  transitions  to  Chapman-Jouguet  (CJ)  detonation  and  the 
wave  velocity  is  constant  thereafter.  The  transition  to  detonation  may  be  smooth  or  it  may  be 
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DISTANCE 


FIGURE  1.  Time-Space  Diagram  of  Runup  of  a  Reactive  Shock  Wave 
to  a  Detonation. 

discontinuous.  The  fraction  of  solid  burned  within  the  reacting  shock  wave  increases  with 
wave  propagation  distance  (or  time)  until  it  is  all  burned  at  Xj.  In  addition  to  combustion, 
which  occurs  instantaneously  within  the  discontinuity,  the  remaining  solid  continues  to  burn 
behind  the  wave  at  some  finite  rate.  The  objective  of  the  flow  field  analysis  is  to  calculate  the 
bum  rate  behind  the  wave  as  a  function  of  pressure  and  the  velocity  of  the  wave. 

To  calculate  the  bum  rate,  one  starts  with  the  Lagrangian  equations  for  conservation  of 
mans,  momentum,  and  energy  (Reference  2) 


*s  *  ~Pm 

(la) 

(lb) 

(lc) 

and  an  arbitrary  equation  of  state 

P  *  P(VJE,W) 

(Id) 

where  V,  U,  E,  P,  and  W  are  the  specific  volume,  velocity,  energy  per  unit  initial  volume, 
pressure,  and  solid  mass  fraction,  respectively.  The  subscripts  t  and  m  denote  partial 
derivatives  with  respect  to  time  and  the  Lagrangian  mass  coordinate.  In  addition,  one  has  the 
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jump  equations  exp  resting  conservation  of  mate,  momentum,  and  energy  across  the  reacting 
shock  wave  (Reference  2) 


pi » 

w  -  W  -  *i> 


£.  a  u\n  (2c) 

and  the  equation  of  state 

Px «  (2d) 

where  the  subscripts  0  and  1  denote  conditions  in  the  uniform  initial  state  and  the  state  just 
downstream  of  the  reactive  shock  wave.  The  velocity  of  the  reactive  shock  is  given  by  (J$. 
Equation  2  assumed  that  Pq  a  o,  Eq  a  0,  and  Uo  a  o.  Note  that  since  the  material  is  allowed 
to  partially  react  within  the  discontinuity,  there  are  six  variables  and  four  equations,  leaving 
two  degrees  of  freedom  instead  of  the  usual  one. 

Forest  calculated  the  burn  rate  Wx  as  a  function  of  thermodynamic  parameters  at  state 
1,  the  pressure  gradient  Pm,  and  what  amounts  to  the  wave  acceleration.  He  found  that  a  zero 
pressure  gradient  Pn  gave  good  results.  For  simplicity,  the  pressure  gradient  is  allowed  to 
vanish  here  early  in  the  derivation. 

The  burn  rate,  Wx,  can  be  obtained  by  differentiating  the  equation  of  state  (Equation  Id) 

P.=V\  +  V!.  +  V\  <3) 

where  the  derivatives 


are  known  functions  of  V,  E,  and  W  and  are  computed  numerically  from  Equation  2d.  The 
time  derivatives  Vx  and  Ex  are  next  eliminated  from  Equation  3  in  favor  of  the  velocity 
gradient,  Um,  using  Equation  1,  yielding 
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*\-C pv-pplum  +  pww< 

It  i«  desired  to  express  the  derivatives  Px  and  UM  in  favor  of  the  pressure  gradient  Pn  and 
derivatives  in  the  direction  of  shock  propagation,  since  these  are  known  at  the  shock 
boundary. 


To  obtain  directional  derivatives,  a  simple  exercise  of  the  chain  rule  gives 


P  +P 

t  m 


(5) 


and 


where  the  notation  (d/ftc),  denotes  a  time  derivative  in  a  direction  on  the  space-time  plane  that 
is  parallel  to  the  direction  of  the  travel  of  the  reactive  shock.  Using  Equation  lb  for  Ux,  while 
noting  that  (dm/3x)t  »  UJV o,  and  evaluating  the  result  at  the  shock  boundary.  Equations  5 
and  6  yield 


P 


i 


-p.u./v* 


(7) 


and 


(8) 


Here  it  is  convenient  to  impose  the  simplifying  assumption  Pm  *  0.  Substituting  Equations  7 
and  8  into  4,  one  obtains 


(9) 


Instead  of  derivatives  with  respect  to  tj,  the  shock  position,  x\,  is  a  more  convenient  variable 
and  so  Equation  9  becomes 


dPx  /  \  dUx 


(10) 


w 


'l  '  — I 

Equation  10  is  a  preliminary  result  obtained  by  Forest  in  Reference  5. 
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In  order  to  relate  the  two  derivatives,  dU\ldx\  and  dP\ldx\,  the  reaction  that  takes  place 
within  the  discontinuity  must  be  specified.  Forest  proposes  a  "reactive”  Hugoniot  in  the  form 
of 


lft  =  Cr  +  Sru  J 

This  equation,  together  with  Equation  2,  determines  at  the  shock  boundary.  Combining 
Equation  11  with  Equation  2a  to  eliminate  the  shock  velocity,  U„  and  differentiating  the 
result,  one  gets 


dPl  /  N  dUt 

oSi.(c+2S^)  — 


Substituting  Equation  12  into  Equation  10,  one  gets 


(12) 


P  —  P  P 
r*  v  1  * 
0  Cr  +  2SrUi 


JLfl 

pw  *1 


relating  the  reaction  rate  IF,  to  the  rate  of  pressure  increase  at  the  shock  front. 


(13) 


An  alternate  and  more  self-consistent  assumption  is  a  nonreactive  shock  front  with 
Wj  =  1.  For  some  equations  of  state,  the  nonreactive  Hugoniot  reduces  to  the  form  of 
Equation  11  with  particular  choices  of  Cr  and  This  is  true  with  the  HOM  equation  of  state 
(Reference  2)  used  by  Forest.  In  this  case,  Equation  13  can  be  applied  to  the  nonreactive  shock 
case  by  proper  choice  of  Cr  and  Sr-  For  other  equations  of  state,  one  can  derive  a  relation 
equivalent  to  Equation  12  as  follows.  First  differentiate  Equation  2d  with  W\  constant  and 
eliminate  the  derivative  dE\!dx\  with  Equation  2c  to  obtain 


dP.  dV.  dU. 


(14) 


Then  combine  Equations  2a  and  2b  to  eliminate  U,  and  differentiate  the  result  to  get 


Eliminate  dV \idx\  between  Equations  15  and  14  and  solve  for 

px-tv*-vi'pv 

*1  =  *  W*  -  V 


(15) 

to  obtain  the  result 

(16) 
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Substituting  Equation  16  into  Equation  10  to  eliminate  dU\!dx\  and  simplifying  the  result 
with  Equation  2b  gives 


3iWy0-VV'V|  <g|  (17) 

Wi  P*  dxl 

which  is  equivalent  to  Equation  13  for  a  npnreactive  shock  but  is  valid  for  any  equation  of 

state. 

Given  the  pressure  rate  of  increase  with  distance,  dPi/dxi,  then  Equation  13  or  17  may 
be  employed  to  calculate  the  required  reaction  rate  (provided  that  Pm  =  0  is  valid).  When  one 
performs  an  experiment  where  Wt  and  the  thermodynamic  state  is  known,  then,  in  effect,  one 
measures  the  reaction  rate  as  a  function  of  that  thermodynamic  state.  Forest  uses  the  wedge 
test  for  this  purpose. 

In  the  wedge  test,  a  shock  wave  is  introduced  into  an  explosive  where  it  accelerates  to 
become  a  detonation  wave  within  a  length  called  the  run  distance.  A  plot  of  the  run  distance 
as  a  function  of  pressure  is  called  a  Pop  plot  after  its  originator,  A.  Popolato.  When  the 
logarithm  of  the  run  distance  is  plotted  against  the  logarithm  of  the  initial  shock  pressure,  a 
straight  line  is  often  obtained  over  a  large  range  of  pressure.  These  data  have  been  obtained 
for  a  wide  variety  of  explosives  and  detonable  propellants.  For  the  wedge  test -experiments 
with  initial  shocks  having  different  initial  pressure -the  "single-curve  buildup  principle”  is 
observed  to  be  a  good  approximation.  (See  Reference  6  for  a  good  discussion  of  this.)  This 
principle  postulates  that  the  shock  that  accelerates  to  detonation  passes  through  a  unique 
space-time  trajectory  regardless  of  the  initial  shock  pressure.  Accordingly,  the  Pop  plot,  taken 
from  the  results  of  several  experiments,  is  also  the  pressure-distance  history  for  any  single 
test  For  a  Pop  plot  correlation  in  the  form  of 

CnXr  =  A  +  B(nP  (18) 

where  A  and  B  are  constants,  Xr  is  the  rim  distance,  and  P  is  the  initial  pressure,  the  single¬ 
curve  buildup  principle  yields 


(19) 


Substituting  Equation  1 )  into  Equation  13  or  Equation  17,  one  gets  an  explicit  expression  for 
the  reaction  rate  as  a  function  of  the  thermodynamic  state  of  the  material. 
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REACTION  KINETICS 


Forest  assumes  that  the  reaction  of  the  explosive  takes  place  according  to  a  simple  first- 
order  decomposition  reaction.  Therefore,  he  fits  the  reaction  rates  calculated  above  to  the 
following  expression  for  first-order  kinetics: 

WW  *  f(P)  (20) 

where  the  reaction  rate  constant  f[P)  is  independent  of  the  solid  mas3  fraction  and  is 
expressed  solely  as  a  function  of  pressure.  The  choice  of  pressure  as  the  independent  variable 
for  the  reaction  rate  constant  is  not  unique;  temperature,  or  internal  energy,  or  any  other 
combination  of  state  variables  (such  as  in  the  H YDROX  computer  code  in  Reference  7)  could 
have  been  just  as  easily  chosen. 

The  first-order  kinetics  formula  is  based  on  the  law  of  mass  action  and  is  useful  as  a 
global  model  of  unimolecular  decomposition  reactions  (Reference  8).  It  is  not  clear  that  this  is 
the  dominating  reaction  mechanism  in  the  shock-initiated  combustion  of  heterogeneous 
explosives.  There  is  a  problem  with  calibrating  Equation  20  at  high  pressures,  where  Wx  is 
finite  and  well-behaved  while  can  become  small  or  vanish,  depending  upon  the  choice  of 
the  reacting  shock  Hugoniot.  In  some  of  the  numerical  hydrodynamics  codes  that  utilize 
Forest  Fire,  the  rate  is  set  to  infinity  whenever  the  pressure  becomes  greater  than  the  CJ 
value.  This  is  accomplished  numerically  simply  by  converting  all  remaining  solid  explosive  to 
gaseous  products  when  the  CJ  pressure  is  exceeded.  This  feature  makes  application  of  Forest 
Fire  difficult  for  problems  where  the  reaction  zone  of  the  detonation  wave  must  be  resolved. 

An  alternate  hypothesis  to  the  first-order  kinetics  assumption  used  by  Forest  is  the 
ignition  and  growth  concept  used  to  describe  explosive  decomposition  by  Lee  and  Tarver 
(Reference  9),  and  Tarver  and  Hallquist  (Reference  10).  The  shock  initiation  of  the 
heterogeneous  solid  explosive  is  modeled  by  ignition  at  localized  hot  spots,  followed  by  grain 
burning  at  the  growing  boundaries.  This  model  has  had  some  success  in  correlating  the 
detailed  mass  velocity  profiles  in  the  flow  field  behind  the  shock  during  shock-to-detonation 
experiments,  lending  credibility  to  the  basic  concept.  According  to  the  ignition  and  growth 
model,  the  global  reaction  rate  is  a  minimum  at  the  shock  boundary  where  W\  ~  1  because  of 
the  jmall  burn-surface  area  of  the  hot  spots.  This  is  in  direct  opposition  to  the  assumption  of 
first-order  kinetics,  where  the  reaction  is  maximum  at  the  shock. 

For  comparison  with  the  first-order  kinetics  assumption,  an  alternate  global  model  for 
the  explosive  decomposition  is 

Wx  =  f{P)  (21) 

which  ignores  all  dependence  of  the  reaction  rate  on  the  solid  mass  fraction.  In  the  following 
calculations,  Equation  21  will  be  referred  to  as  the  constant  rate  model.  This  model  should  lie 
somewhere  between  the  extremes  of  the  ignition  and  growth  burn  model  and  the  first-order 
kinetics  employed  by  Forest.  It  can  be  noted  that  Equation  21  does  not  have  the  singularity 
when  VPj  =  0,  as  Equation  20  does. 
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NUMERICAL  CALCULATIONS 


The  effect  of  the  various  assumptions  employed  in  the  Forest  Fire  burn  model  derivation 
are  best  tested  numerically.  The  calculations  will  employ  the  explosive  PBX-9404  because  of 
the  large  amount  of  experimental  data  available  for  it.  The  equations  and  constants  for  the 
HOM  equation  of  state  for  PBX-9404  are  given  in  Reference  2.  These  are  used  except  for  the 
Hugoniot  of  the  unreacted  explosive  and  the  wedge  test  run  distance  correlation  (Pop  plot), 
which  are  obtained  from  Reference  1 1. 

Partially  reacting  Hugoniot  curves  with  constant  l?i  are  calculated  for  9404  and  ore 
plotted  on  the  P-V  plane  in  Figura  2.  The  unreacting  shock  Hugoniot  is  the  curve  labeled 
ffj  a  1  in  the  figure,  and  the  CJ  point  for  steady  detonation  is  labeled  on  the  combustion 
products  (W\  =  0)  curve.  Reacting  Hugoniot  curves  corresponding  to  different  choices  of  the 
constants  Cr  and  Sr  in  Equation  11  can  also  be  represented  on  the  P-V  plane.  For  one  choice  of 
the  constants,  the  reacting  Hugoniot  passes  through  both  the  initial  state  and  the  CJ  point. 
This  Hugoniot  is  shown  on  Figure  2  and  will  be  used  for  the  reacting  shock  example  in  the 
following,  calculations.  The  explosive  solid  fraction  Wi  is  not  constant  on  this  curve.  The 
"ariation  of  solid  fraction  with  pressure  along  the  reacting  shock  Hugoniot  is  shown  in  Figure 
3.  The  unreacting  shock  Hugoniot  that  will  be  used  for  comparison  in  the  calculations  is  the 
Wy  a  1  curve  *j i  Ficure  2. 

The  reaction  rcte,  was  calculated  using  Equation  13  for  both  reacting  and 
unreacting  shock  Hypotheses,  and  the  result  plotted  in  Figure  4.  At  low  pressures,  the  results 
are  nearly  identical.  It  it  only  at  pressures  near  the  CJ  point  that  the  reactive  shock  case 
gives  slightly  higher  reaction  rates.  The  difference  between  the  two  curves  is  negligible 
compered  to  differences  induced  by  the  choice  of  the  reaction  model. 

In  Figure  5.  the  quotient  to  plotted  as  a  function  of  pressure  for  both  the  reacting 
shock  case  arid  the  nonreacting  case,  where  for  the  letter  Wi  =  1.  The  difference  between  the 
two  results  is  suostaiitial,  partto  -Jarl*'  at  high  pressures  where,  according  to  Figure  3, 
goes  to  xerc  ns  the  CJ  pressure  is  approached.  Physically,  very  little  reaction  takes  place 
within  the  shock,  and  the  numerical  codes  where  the  Forest  Fire  burn  model  is  employed  do 
not  include  such  a  phenomenon.  Therefore,  the  assumption  of  a  reacting  shock  wave  used  by 
Forest  to  calculate  the  reaction  rates  is  not  satisfied  in  the  resulting  numerical  calculations. 
The  further  .wumpiio.i  of  first-order  kineHca  exaggerates  the  consequences  of  the  reactive 
shock  assumption. 

To  compare  the  effect  of  the  various  assumptions  and  models,  calculations  are  made  of 
the  run  distance  obtained  in  the  wedge  test.  A  one-dimensional  Lagrangian  hydrodynamics 
computer  code  was  utsed  for  the  calculations.  It  uses  the  SIN  code  algorithm  described  in 
Reference  2. 

There  is  some  arbitrariness  in  a  general  simulation  of  the  wedge  test  because  of  the 
different  types  of  shock  attenuators  and  explor;  ve  drivers  actually  used  in  the  tests.  These  are 
necessary  in  order  to  introduce  shock  waves  of  different  intensities  into  the  explosive 


10 


NWC  TP  6898 


J 


jDq>j  ‘ajnssajj 

fl 


n 


1 


KIGURK  2.  Partially  Reacting  ( If  j  =  Constant)  Hugnniol  ('urves  for  PBX-9404  The  routine 
shock  I  lugonict  that  passes  through  the  initial  and  CJ  states  is  shown. 


NWCTP  6898 


uojpojj  pj|os 


12 


FIGURE  3.  Solid  Fraction  as  a  Functioning  of  Shock  Pressure  for  Reactive  llugoniot  Plotted 
in  Figure  2. 
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PIGURK  4.  Reaction  Rate  -  IV,  Plotted  as  Function  of  Pressure  Comparing  Results  (Mourned 
With  Reactive  llugoniotund  Nonreactive  Shock  Assumptions. 
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FIGURE  5.  Reaction  Rate  -  H\/W|  Showing  Effect  of  First-Order  Kinetics  Assumpti 
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material.  Calculation*  show  that  tha  run  distenc*  in  tha  explosive  depends  on  the  shock 
imped*"*?*  of  tha  attenuators.  To  evaluate  tha  magnitude  of  this  affect,  Figure  6  shows  run 
dlftimni  calculated  as  a  function  of  initial  shock  pressure  for  two  extreme  cases.  In  the  first 
case,  the  shock  was  produced  by  the  impact  of  a  plexiglas  flyer  plate.  In  the  second,  the  shock 
was  generated  by  imposing  a  constant  velocity  (a  piston)  upon  the  boundary  of  the  explosive 
material,  thus  simulating  an  infinite  shock  impedance.  The  HOM  equation  of  state  for 
pfoidgi***  (Reference  2)  was  used.  For  this  calculation,  constant  rate  kinetics  (Equation  21) 
wee  used  for  the  burn  model  and  was  fitted  using  the  nonreacting  shock  Hugoniot.  The  run 
distance  for  the  plexiglass  flyer  plate  in  Figure  6  is  approximately  20%  longer  than  the 
corresponding  distance  calculated  for  the  constant  velocity  piston.  In  the  following 
calculations,  the  plexiglass  flyer  is  used. 

The  influence  of  the  reaction  model  on  run  distance  is  shown  in  Figure  7.  The  run 
^i«»anrw  normalized  by  the  experimental  value  obtained  from  the  Pop  plot  is  plotted  against 
the  input  shock  pressure.  Each  of  the  two  kinetics  models  is  fitted  to  the  extreme  choices  of 
nonreacting  shock  Hugoniots  and  the  reacting  Hugoniot  which  passes  through  the  CJ  point. 
The  curves  are  labeled  in  the  figure  for  first-order  kinetics: 

FFl-reacting  Hugoniot 
FF2-nonreacting  shock 

and  for  constant  rate  kinetics: 

FF3-reacting  Hugoniot 
FF4-nonreacting  shock 

For  the  first-order  kinetics  model,  a  wide  range  of  response  can  be  obtained  by  the  choice  of 
the  reacting  Hugoniot  The  calculated  results  at  low  initial  shock  pressure,  however,  are  not 
very  sensitive  to  the  kinetic  model  or  its  calibration.  This  is  because  the  solid  mass  fraction  is 
nearly  unity  at  low  pressures  on  the  reacting  Hugoniot  used  for  the  calibration,  and  this 
condition  is  obtained  for  most  of  the  run  distance. 

The  difference  in  behavior  of  the  different  models  is  most  extreme  at  the  high-pressure 
end.  This  pressure  region  in  of  importance  in  physical  problems  where  the  detonation  reaction 
zone  must  be  resolved.  These  problems  include,  for  example,  critical  diameter,  response  to 
thin,  high-speed  flyer  plates,  and  minimum  priming  charge  tests.  Any  calculation  of  these 
types  of  problems  using  a  variation  of  the  Forest  Fire  burn  model  should  be  viewed 
skeptically  until  thorough  experimental  verification  is  obtained. 

On*  consequence  of  the  constant  rate  kinetics  model  is  that  its  implementation  is  more 
automatic  than  the  first-order  kinetics  model.  The  range  of  calculated  run  distance  shown  in 
Figure  7  is  much  less  between  tha  r« acting  and  nonreacting  choices  of  Hugoniot  This  means 
that  one  can  get  a  good  correlation  with  the  wedge  test  experiment  with  a  minimum  of 
guessing  at  reacting  Hugoniots.  In  this  regard,  however,  it  can  be  noted  from  Figure  7  that 
the  constant  rate  curves,  FF3  and  FF4,  are  reasonably  straight  lines,  so  some  improvement 
could  be  obtained  by  artificially  modifying  the  input  Pop  plot. 
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PRESSURE,  kbar 

FIGURE  fi  Wedge  Test  Simulation  Showing  Effect  of  Flyer  Plate  Impedance. 
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FIGURE  7,  Kll't  vl  Reactive  Hugoniot  Assumption  and  Reaction  Kinetics 
Model  on  Wedge  Test  Simulated  by  a  Plexiglas  Flyer  Plate. 
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In  the  course  of  making  the  calculations,  plots  of  pressure  profiles  of  the  reactive  flow 
were  made.  Initially  the  pressure  gradient  was  quite  small  and  Pn  =  0  appeared  to  be  a  valid 
assumption.  However,  at  later  times,  after  a  period  of  shock  acceleration,  a  very  strong 
gradient  developed.  This  being  the  case,  it  is  not  clear  why  the  Forest  Fire  burn  model  gives 
such  good  overall  agreement  for  the  wedge  test  run  distance.  If  one  examines  individual 
experimental  wave  trajectories  corresponding  to  different  initial  shock  pressure  to  verify  the 
single-curve  buildup  principle,  one  finds  curves  of  wave  velocity,  U9,  plotted  against  wave 
position,  X\,  which  all  lie  in  a  narrow  band.  However,  individual  shock- wave  accelerations 
within  this  band  can  deviate  substantially  from  the  mean  and  from  one  another.  One  can 
speculate  that  input  of  an  average  shock  trajectory  will  yield  a  good  prediction  of  global 
results  like  run  distance,  even  though  details  of  the  process  may  be  erroneous. 


CONCLUSIONS 


The  reaction  rate,  Wt,  expressed  as  a  function  of  pressure,  is  not  sensitive  to  the  choice  of 
reactive  or  nonreactive  Hugoniot  used  in  its  calibration.  The  choice  becomes  important  when 
the  formula  for  first-order  kinetics  is  assumed  to  express  the  rate.  With  the  first-order 
kinetics  formula  calibrated  first  using  the  extreme  reactive  and  then  the  nonreactive 
Hugoniot,  the  calculated  wedge  test  results  bracket  the  experimental  results.  One  can 
therefore  fine-tune  the  agreement  with  experiment  by  adjusting  the  choice  of  reactive 
Hugoniot  The  drawback  to  this  approach  is  that  neither  the  first-order  kinetics  nor  the 
reactive  Hugoniot  has  any  physical  relation  to  shock-initiated  burning  in  heterogeneous 
explosive.  A  simpler  constant  reaction  rate  approach  is  suggested  that  avoids  the  difficulties 
of  choosing  an  unphysical  reactive  Hugoniot  while  giving  a  satisfactory  correlation  with  the 
wedge  test  results. 

Calculations  of  wedge  test  results  at  low  initial  shock  pressure  were  basically 
insensitive  to  the  kinetic  model  and  the  Hugoniot  used  in  its  calibration.  The  differences  were 
quite  large,  however,  for  wedge  test  calculations  starting  with  high  initial  shock  pressures. 
The  reason  is  that  at  low  pressure,  the  solid  mass  fraction  is  nearly  one  for  most  of  the  run 
distance.  Under  this  circumstance,  there  is  no  essential  difference  between  the  various 
reaction  models  and  their  calibration. 

The  Forest  Fire  burn  law  as  formulated  by  Forest  is  not  applicable  to  high-pressure 
phenomena  where  the  pressure  may  exceed  the  CJ  value,  this  situation  is  obtained,  for 
example,  whenever  the  reaction  zone  in  a  detonation  must  be  resolved.  The  constant  rate 
approach  discussed  here  may  be  a  reasonable  way  to  extrapolate  moderate  pressure  reaction 
rates  to  pressures  greater  than  the  CJ  value.  The  validity  of  this  approach  needs  to  be  tested. 
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